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Abstract 

In relativistic heavy-ion collisions, the system has gone through a series of evolution, almost 
at every stage of its evolution it leaves behind footprints in flow observable. Those footprints 
contain valuable information of the bulk property of the (nearly) perfect liquid. By examing 
footprints of the nearly perfect liquid, we address a few important issues, including the ideal 
hydrodynamic limit, estimation of 77/s, testing the Number of Constituent Quark scaling at low 
energy, in small system, at large transverse momentum, and in forward region. Future prospect 
of flow study is discussed. 



1. Introduction: the (nearly) perfect liquid 

As a unique tool to study the QCD matter under extreme conditions, RHIC has been success- 
ful in operations since year 2000. The wealth of data, collected and analyzed in many aspects, 
indicates that central Au+Au collisions can be well described by ideal Hydrodynamics yj]. In- 
dications of liquid-like behavior of the matter that RHIC has created came in the form of large 
elliptic flow. What is more interesting is that, this liquid has little viscosity and acts like a perfect 
one yfl. Those findings lead to the announcement of the discovery of the existence of a perfect 
liquid!!. 

In following sections, we address a few important questions which are directly related to the 
understanding of the property of the (nearly) perfect liquid. 

2. Is hydrodynamic limit saturated? 

The search for fluid-type signature started 30 years ago at BEVALC. Till 2000, when RHIC 
started operation, data becomes close to hydrodynamic predictions. This can be seen by plotting 
V2 scaled by the initial eccentricity e as a function of 1/S dN/dy, data points reach hydro limit 
in most central collisions fl. However, as our understanding advances, it is realized that when 
calculating the hydrodynamic limit, there is ambiguity in the choice of initial condition @|, 
uncertainties in the Equation of State (EoS) [6], and hadronic dissipative effects are also not 
negligible [Q]. Thus it is important to revisit the issue. 

In Figure[TJ V2{4), scaled by the standard eccentricity e st d, is plotted as a function of transvers 
momentum (p,) for difference centralities. We see that vz/s increases with centrality, due to the 
increase of 1/S dN/dy. For V2 at a fixed p t , one finds that V2 keeps increasing until it saturates at 
the bound set by hydrodynamics. This explains the saturation of p t integrated V2 in Jit], without 
the need for a softness in the EoS. 1 . 



1 As a related remark, Heinz and Kestin pointed out at this meeting that the non-monotonic feature of V2 at a fixed p, 
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Figure 1: V2 measured by four particle cumulant (i'2{4)) scaled by standard eccentricity used in hydro calculations, as 
a function of transverse momentum (p t ). Data points are for Au+Au collisions at 200 GeV. Curves are hydro dynamic 
predictions. Data points of this plot is from @|. 
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Figure 2: 1'2/s scaled by the corresponding saturation value (hydro limits) obtained from the simultaneous fitting, for 
Glauber (left) and CGC (right) initial conditions. The hydro limit is by definition centered at unity, with error represented 
by the shaded bars. Data points are for Au+Au collisions at 200 GeV. 
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Figure 3: r//s as a function of 1 /S dN/dy for collisions at 200 GeV. The conjectured quantum limit, as well as rj/s for 
He at T c is also plotted for comparison. This plot is from fl4ll 



For the p, integrated v%, to quantify the possible discrepancy between data and ideal hydro- 
dynamics, one can either match the data using hydrodynamic models which incorporate viscous 
U], or fit the data with a formula motivated by transport models 1 12, 131 



corrections 

Figure |2 shows V2/S scaled by its corresponding limit at which they saturates (hydro limit). 
The hydro limit is obtained from fitting vi/e with the transported model motivated formula 
— = — , , J,„ , where K is the Knudsen number defined by the mean free path A scaled by 

£ L £ Jhydro 1+K/Ko' ■> f J 

the system size R, and Kq is obtained from transport calculations. This formula has the desired 
feature of describing the system at two extremes; when K is small, the deviation to the saturation 
value (hydro limit) is proportional to K, which corresponds to first order correction of viscosity, 
and when K is large, proportional to 1/K, which corresponds to low density limit. A mass hi- 
erarchy is observed, the heavier the particle, the more saturation is seen. For charged particles, 
we see that the system is still 30-50% away from ideal hydro dynamics, depending on initial 
conditions. 

The Knudsen number obtained from this procedure can be used to calculate 77/s. Follow- 
ing Il5ll . the viscosity for a classical gas of massless particles with isotropic differential cross 
sections is 77 = 1.264T/<x llrjll . It is arguable to apply the formula to strongly interacting 
dense matter, however, in practice the viscosity recovered from this procedure agrees well with 
that obtained from viscous hydro calculations 11411 . Taking the entropy density for a classi- 
cal ultra-relativistic gas as s = 4n, with n the particle density, then 77/s can be calculated as 
77/s = 0.316KRT. The temperature T is obtained from fitting STAR'S pion mx slope Intl . In 
Figure [3] 77/s is plotted as a function of 1/S dN/dy for Glauber and CGC initial conditions. For 
both 77/s is lower than that for He at T c . 77/s for CGC initial condition is smaller than that for 
Glauber initial condition, because with CGC initial condition, a stronger saturation is seen in 
the shape of vt/s vs. 1/S dN/dy, which gives a smaller K for same 1/S dN/dy. This does not 
necessarily contradict to the conclusion arrived from viscous hydro calculations J^lQj]], m which 
the Equation of State is chosen to be the same for the two initial conditions. 



, as shown in floll . is related to the interplay between radial flow and freeze-out. It can not be associated unambiguously 
with a phase transition in the EoS. 

3 



a>.E 
o 5 

03 B 



(CGC) . 



(Glauber) 



PHENIX v 4 WWND09 
Lacey et al, PRL 092301 (2007) 



A * Drescher et al. PRC76 024905 (2007) 
(CGC) (Glauber) 



STAR p t corrosion (2009) 



Gavin & Abdel-Aziz. PRL 97 162302 (2006) 

(p t correlation) £ (number density correlation) 

Hydro, calculations, Song SQM08, Heinz WWND09 
P. Romatschke & U. Romatschke, PRL 99 172301 (2007) 
PHEN IX PRL 98 172301 (2007) 

Hees et al., arXiv:0808.3710 



KMeyer, PRD 76, 101701(R) (2007) [Lattice QCD] 
^■^P^b> (T = 1 .65 T ) 

Demir & Bass, arXiv:081 2.2422 ( 2009) [had ron gas] 

__l I i i i I i i i I i i 



STAR Chg. V (2009) 



8 



10 12 

47t T|/S 



Figure 4: Collection of T]/s calculations. Color code : Red for results based on flow measurements, cyan for results 
based on fluctuation measurements, orange for results based on hydrodynamic calculations, blue for estimation based on 
heavy flavor measurements, and magenta for other calculations. 



3. What is the perfection ? 

The 77/s presented in previous section is an effective quantity which includes viscous effects 
over different phases. It is desirable to extract 77/ s in QGP fluid phase, however, all "signatures", 
spectra, flow, fluctuations, involve some time average over the history of plasma. One has to rely 
on models to infer the 77/s from QGP phase. It thus becomes important to understand the relative 
contribution from phases other than the QGP fluid phase. Theoretical work llTci.[l9ll indicate that 
above T c , the 77/s increases with T, and below T c , decreases with it (Fig. [5]). The divergence 
happens at T c and one can infer a minimum of 77/s near it. Note that if the system spends even 
a short period of time in the hadronic phase it may catch considerable amount of viscous effect. 
Indeed it is calculated lEoll that the 77/s for hadronic gas is 0.4-0.5, a few times of the conjectured 
quantum limit. Nevertheless, if one collects 77/s calculations on the market (Fig. 3J, although the 
result spread in a wide range from to 10 times of the conjectured quantum limit, which reflects 
the current uncertainty in the calculation of 77/s, most of the calculated 77/s is below that of He at 
T 



4. What is the ultimate say on partonic collectivity ? 

<p and Q have been suggested as an ideal probe for partonic collectivity, because they have 
small hadronic cross section thus carry information from the early, partonic stage. The large 
statistics accumulated during RHIC run VII makes it possible to measure V2 for </> and Q. with 
unprecedented accuracy. In the right panel of Figure |6l it is found that </> and Q. has sizable vi 
and they still follow Number of Constituent Quark (NCQ) scaling, just like light hadrons (shown 
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Figure 5: (left) 77/s calculated for low energy pion gas and quark-gluon plasma fT3l . (right) tj/s calculated from weakly 
coupled real scalar field theories fish . 
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Figure 6: v% fo r lig ht hadrons (left) and multi strange hadrons (right). Open symbols is from Phenix [2 ill and solid 
symbols, STAR [22]. The curves are NCQ motivated fits (23ll . 



in the left panel for comparison). This is the definite proof that partonic collectivity has been 
reached at RHIC. 



5. What is the limit of Number of Constituent Quark scaling ? 

The flow pattern of baryons and mesons can be explained well by the Number of Quark 
scaling [24], which is viewed as an evidence for the existence partonic degree of freedom. One 
does not expect the partonic phase to exist when the system is small and/or when the energy of 
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the system is low, thus as a consequence, the NCQ scaling is expected to break down eventually. 
As a controlled reference, it would be interesting to find out the condition under which the NCQ 
scaling breaks down. 
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Figure 7: v% scaled by number of quarks (n q ) as a function of kEj /n„, for Pb+Au collisions at 158 AGeV 12511 (left) and 
CuCu collisions at 200 GeV (right). 



In the left panel of Figure [7] V2 of identified particles are scaled by their constituent quarks 
and plotted as a function of kEjln q , for Pb+Au collisions at 158 AGeV. We see that within 
errors, the result is still consistent with NCQ scaling up to kEjln q = 0.8GeV/c 2 . In the right 
panel, a similar plot is shown for Cu+Cu collisions at 200 GeV, V2/n q for identified particles fall 
in a common trend well. No sign shows that NCQ breaks at these two conditions. NCQ is also 
studied at forward region in M26I1 . within error there is no sign of breakdown. 

Besides the partonic degree of freedom, NCQ scaling also implies that particles at interme- 
diate p, are produced by quark coalescence. At large p r as hard process begins to kick in and 
particles are no longer produced by quark coalescence, thus the NCQ scaling is expected to break 
down. It is important to locate the p, range where the NCQ scaling breaks down, as it tells us 
the region of transition of particle production mechanism. It is as well important to examine 
the pattern with which the NCQ scaling breaks down for various hadrons. Such pattern will not 
only shed a light on the dynamics of jet fragmentation, it will also, being served as a counter 
example of NCQ scaling, deepen our understanding of quark coalescence. Figure[8]shows V2/n q 
for identified particles from both PHENIX (left) and STAR (right). It is obvious that \'2jn q of 
protons and pions begin to diverge starting around kEx/n q = lGeV, possible difference between 
lambda and pions is also observed. 



6. Outlooks 



From the (3+l)-dimensional hydrodynamics with exact longitudinal boost-invariance 1280, 
it is demonstrated that a phase transition from hadronic matter to QGP leads to non-monotonic 
behavior in both beam energy and impact parameter dependence. This is shown in right panel of 
Figure |9j in which for a combined Equation of State (EOS Q) that includes the phase transition 
between a resonance gas (EOS H) and for the QGP phase (EOS I), a dip is observed at the range 
of SPS energy. As vo begins to rise again for high density at LHC energies, the dip only covers 
the energy range between SPS and RHIC. The significance and the range it happens has to be 
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Figure 8: V2 scaled by number of quarks (n q ) as a function of kEjfn q [27|](left) and p t /n q 12211 (right), for Au+Au 
collisions at 200 GeV 



determined by experiments. Left panel of Figure [9] shows the excitation function of V2. We see 
a general trend of v% increase with energy (thus density), with current error we cannot identify 
any possible abnormal around -y/iNN=10 Gev, the same energy at which the k + /n + "horn" 12911 
is located. It is desirable to return to this energy range with accurate measurement to identify 
possible abnormalities. 

Looking into RHIC II program and future heavy ion program at LHC, as part of a compre- 
hensive test for the deviation from hydro limit, it is important to i) increase the system size while 
keep the energy the same (U+U collisions), and ii) increase the energy while keep the system 
size the same. The former will allow us to test if V2 saturates or not by increasing the particle 
density in the transverse plane, and the latter will allow us to check if the viscosity increases as 
anticipated once we passed the transition temperature. 
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Figure 9: Left: V2 excitation function, the plot is from 1 30]. Right: Hydro calculation of integrated V2 as a function of 
pion multiplicity, for a given impact parameter. The plot is from |28|| 



7. Summary 

The deviation from ideal hydrodynamic limit is discussed. From the transport model based 
approach, it is found that the system is 30-50% away from the value at which it is supposed to 
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saturate. The perfection of the nearly perfect liquid is reviewed. The current estimations of 77/s 
spread in a wide range from to 10 times of the conjectured quantum limit, however, most of 
the calculated 77/s is below that of He at T c . It is concluded that the partonic collectivity has been 
reached at RHIC. The number of quark scaling is found to hold roughly at low energies, in small 
systems and in forward regions, however, signs of break down of NCQ scaling are observed 
at large p,. To the aspect of understanding the bulk property of the nearly perfect liquid, the 
importance of future heavy ion programs at RHIC and LHC is emphasized. 
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